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Abstract: 

Understanding the formation and evolution of the first stars and galaxies represents one of the most 
exciting frontiers in astronomy. Since the universe was filled with neutral hydrogen at early times, the 
most promising method for observing the epoch of the first stars is using the prominent 2f-cm spectral 
line of the hydrogen atom. Current observational efforts are focused on the reionization era (cosmic 
age t ~ 500 Myr), with earlier times considered much more challenging. However, the next frontier 
of even earlier galaxy formation (t ~ 200 Myr) is emerging as a promising observational target. This 
is made possible by a recently noticed effect of a significant relative velocity between the baryons and 
dark matter at early times. The velocity difference suppresses star formation, causing a unique form 
of early luminosity bias. The spatial variation of this suppression enhances large-scale clustering and 
produces a prominent cosmic web on f00 comoving Mpc scales in the 21-cm intensity distribution. This 
structure makes it much more feasible for radio astronomers to detect these early stars, and should 
drive a new focus on this era, which is rich with little-explored astrophysics. 

Keywords: Write keywords here galaxies: formation — galaxies: high-redshift — intergalactic 
medium — large-scale structure of universe — cosmology: theory 



1 Introduction 

1.1 Cosmic reionization and the 
first stars 

Galaxies around us have been mapped system- 
atically out to a redshift z ~ 0. 3 by recent 
large surveys (|Aihara et alj 120111 : IColless et al.l 
2001). The observed galaxy distribution shows 
a large-scale filament-dominated "cosmic web" 
pattern that is re produced by cosmological numer - 
ical simulations (|Springel. Frenk. fc White! [2006). 
This structure is well-u nderstood theoretically 
(|Bond. Kofman. fc PogosvarJ 1 19961 ) as arising from 
the distribution of the primordial density fluctuations, 
which drove hierarchical structure formation in the 
early universe. Recent observations have been pushing 
a new frontier of early cosmic epochs, with individual 
bright gala xies detec t ed re liably from as early as 
z — 7.2 l|Ono et al.1 I2012I ). which corresponds to 
t ~ 750 Myr after the Big Bang. 

The refurbished Hubble Space Telescope [HST] 
and the Spitzer Space Telescope have together pro- 
vided increasingly detailed spectral information on 
populations of z ~ 6 galaxies, making possible robust 
determinations of their stellar masses, though substan- 
tial uncertainties remain in the derived s tellar ages 
and, thus, their star- forma t ion h i stories (|Ono et al 



2010; ISchaerer fc de Barrosl l201Ct ICurtis-Lake et al 



2013). The James Webb Space Telescope [JWST], 



planned for launch in 2018, will have a 6.5 meter mir- 
ror a nd focus on infrared wavelengths l|Gardner et al.l 
2006). This should allow it to reach much higher 
redshifts, detecting dwarf galaxies and perhaps super- 



novae or gamma-ray burst afterglows. 

A major goal of current and upcoming observations 
is to probe the era of cosmic reionization. Ever since 
the discovery that the inte rgalactic gas throughou t the 
Universe is highly ionized l)Gunn fc Petersonl 19651 ) , as- 
tronomers have been searching for the period when the 
hydrogen was ionized for a second time after it became 
neutral at cosmic recombination. This reionization 
is the most recent cosmic phase transition, whereby 
the gas was heated and ionized throughout the uni- 
verse, affecting subsequent galaxy formation and po- 
tentially detectable through its effect on a large range 
of observations. While the highest-redshift galaxies 
may be probing the late stages of reionization, there is 
also a global constraint from the re-scattering of some 
of the photons of the cosmic microwave background 
(CMB) by the freshly created plasma. According to 
the WMAP and Planck CMB satellit es, the bulk of 
cosmi c reionization occurred at z ~ 11 IjHinshaw et al] 
l2012l : lPlanck Collaboration et alj|20l3 ). 

While this observational progress has been remark- 
able, theory suggests that a rich, varied history of 
early galaxy formation still remains to be discovered. 
In particular, current and planned telescopes can di- 
rectly detect early galaxies only if they are particu- 
larly massive or contain a bright, rare event. How- 
ever, it is theoretically expected that the bulk of the 
early stars formed in a large number of very small 
galactic units, which would be difficult to observe indi- 
vidually. Even the requirement of cosmic reionization 
likely implies the exist ence of smaller ga l axies than cur- 
rently observed [e.g., iRobertson et al.1 l|2013T )]. Going 
to even earlier times, high-resolution numerical simula- 
tions suggest that the truly earliest stars formed within 
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~ 10 6 Mp) dark matter halos J A bel, Br yan, fe Normanl 
2002; iBromm, Coppie. fc Larson! 119991 . Because this 
characteristic mass of early galactic host halos is so 
small, constraining their abundance observationally 
would probe primordial density fluctuations on ~ 20 
kpc scales, an order of magnitude below current con- 
straints. This could lead to new limits on models with 
suppressed small-scale power such as wa rm dark mat- 
ter (JBarkana. Haiman. fc Ostr ikcr 2001). 

The first-star simulations mentioned above can 
only follow small cosmic volumes, and thus begin 
to form stars much later than in the real universe 
(since the simulations do not sample the over-dense 
regions that are the sites of the earliest star forma- 
tion), but analytical methods show that the very first 
stars must have formed at z ~ 65 (age t ~ 35 Myr ) 
(jNaoz. Noter. fc Barkana! |2006| ; iFialkov et all I2012T ). 
The formation of the very first star was a major mile- 
stone which ended the dark ages of the universe - the 
period after cosmic recombination, when the universe 
was dark except for the fading glow of the CMB. It also 
marked a transition point from the simplicity of the 
early universe, which was homogeneous and isotropic 
except for small fluctuations that can be described by 
linear perturbation theory. The formation of the first 
stars initiated a new era of complexity and feedback, 
whereby the nuclear and radiative processes within 
stars (on the scale of 10 11 cm) began to affect the 
global state of the cosmic gas on scales larger by up to 
15 orders of magnitude. 

1.2 21-cm cosmology 

The best hope of observing the bulk population of early 
stars is via the cosmic radiation fields that they pro- 
duced. The mean radiation level traces the cosmic 
star formation rate, while spatial fluctuations reflect 
the clustering of the underlying sources, and thus the 
mass of their host halos. In particular, the hyperfine 
spin-flip transition of neutral hydrogen (H I) at a wave- 
length of 21 cm is potentially the most promising probe 
of the gas and stars at early times. Observations of this 
line at a wavelength of 21 x (1 + z) cm can be used to 
slice the universe as a function of redshift z and obtain 
a three-di mensional map of the diffuse H I distribution 
within it fHog an fc Reesl [19791 ) . in the previously un- 
explored era of redshifts ~ 7 — 200. 

Absorption or emission by the gas along a given 
line of sight changes the 21-cm brightness temper- 
ature Tj,, measured relative to the temperature of 
the background source, whic h here is the CMB 
(|Madau. Meiksin, fc Ree3ll997l ). The observed T b is 
determined by the spin temperature T s , an effective 
temperature that describes the relative abundance of 
hydrogen atoms in the excited hyperfine level com- 
pared to the ground state. Primordial density inho- 
mogeneities imprinted a three-dimensional power spec- 
trum of 21-cm intensity fluctuations on scales down 
to ~ 10 kpc (all sizes co moving), making it the rich - 
est data set on the sky (|Loeb fc Zaldarriagal 12004 ). 
In particular, since 21-cm measurements can in prin- 
ciple reach such small scales and do it over a sub- 
stantial redshift range, they can access more inde- 



pendent modes, e.g., a factor of ~ 10 9 more than 
the CMB, making 21-cm measurements potentially 
much more sensitive to small levels of non-Gaussianity. 
The potential yield of 21-cm observations is further 
enhanced by th e expected an isotropy of the 21-cm 
power spectrum (Barkana & Locb 2005a; Nusscr 2005; 
lAli. Bharadwai. fc Pandevl l2005; Barkana 2006a|). 

The 21-cm signal vanished at redshifts above z ~ 
200, when the gas kinetic temperature, Tk, was close 
to the CMB temperature, T 7 . Subsequently, the gas 
cooled adiabatically, faster than the CMB, and atomic 
collisions kept the spin temperature T s of the hyperfine 
level population b elow T-,, so that the gas appeared in 
21-cm absorption l|Scott fc Reesll 19901 1. As the Hubble 
expansion continued to rarefy the gas, radiative cou- 
pling of T s to T 7 started to dominate over collisional 
coupling of T s to Tk and the 21-cm signal began to 
diminish. 

Once stars began to form, their radiation produced 
feedback on the intergalactic medium (IGM) and on 
other newly-forming stars, and substantially affected 
the 21-cm radiation. The first feedback came from the 
ultraviolet (UV) photons produced by stars between 
the Lya and Lyman limit wavelengths. These photons 
propagated freely through the universe, redshifted or 
scattered into the Lyq resonance, and coup led T s to 
Tk once again ()Madau. Meiksin. fc Reesll997[) through 
the Wouthuvsen-Field (|Wouthuvsenlll952l ; lFieldlll958f ) 
effect by which the two hyperfine states are mixed 
through the absorption and re-emission of a Lya pho- 
ton. Meanwhile, Lyman- Werner (LW) photons dissoci- 
ated molecular hydrogen and eventually ended the era 
of p rimordial star forma tion driven by molecular cool- 
ing (JHaiman et al.lfl997D . leading to the dominance of 
larger halos. X-ray photons also propagated far from 
the emitting sources and began early on to heat the gas 
(JMadau. Meiksin. fc Reeslll997l ). Once T s grew larger 
than T 7 , the gas appeared in 21-cm emission over the 
CMB level. Emission of UV photons above the Ly- 
man limit by the same galaxies initiated the process of 
cosmic reionization by creating ionized bubbles in the 
neutral gas around these galaxies. See Figure [T] for a 
brief summary of early cosmic history. 

Several arrays of low-frequency radio telescopes 
are currently being constructed in order to detect the 
21-cm fluctuations from cosmic reionization. Current 
efforts include the Murchison Wide- field Array [MWA 
(JBowman. Morales, fc Hewittl l2009f)l. the Low Fre- 
auencv Array [LOFAR l|Harker et al.ll2010])]. the Giant 
Metrewave Radio Telescope [GMRT (jPaciga et all 
1201 if )], and the Precision Array to Probe the Epoch 
of Reionization [PAPER (|Parsons et al.ll2010D ], and 
early plans have been m ade for a future Square 
Kilom eter Array [SKA; e.g.. lFurlanetto. Oh. fc Briggsl 
(2006)]. Although the expected foregrounds (dom- 
inated by Galactic synchrotron) are much brighter 
than the 21-cm signal, they are not expected to 
include sharp spectral features. Thus, the prospects 
for extraction of the 21-cm signal (and from it the 
reionization history) ar e quite pr omising, using the 
21-cm power spectrum dBowman. Morales, fc Hewittl 
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Figure 1: Overview of cosmic history, with the age of the universe shown on the top axis and the corre- 
sponding redshift on the bottom axis. Yellow represents ionized hydrogen, and gray is neutral. Observers 
probe the cosmic gas using the absorption of background light (dotted lines) by atomic hydrogen. Stars 
formed in halos whose typical size continually grew with time, going from the first generation that formed 
through molecular-hydrogen cooling (red dots) , to the larger galaxies that formed through atomic cooling 
and likely dominated cosmic reionization (blue dots), all the way to gala xies as massive as the Milky Way, 
some of which host bright quasars (green dots). From LBarkanal ( 2006bl ). 
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20051: iFurlanetto. Zaldarriaga. fc Hernauisd l2004bl : 
Ichikawa et alj|20ld : iPan fc Barkanall2012T l. A differ- 



ent approach is to measure the total sky spectrum 
and detect the global reionization signal arising 
from the overa ll disappeara nce of atomic hydrogen 
(iBowman fc Rodeerd l2010l : iPritchard fc Loebl l20ld : 
iMorandi fc Barkanall2012 



2 Large spatial fluctuations 
in galaxy numbers at high 
redshift 

A broad, common thread runs through much of the 
theoretical development of early galaxy formation over 
the last decade: The density of galaxies (or stars) 
varies spatially, with the fluctuations becoming sur- 
prisingly large at hi gh redshift, even on qu ite large 
cosmological scales ( Barkana fc Loebl 120041 '). This 
can be unders tood from the standa r d theory of 
galaxy biasing dPress fc Schechterlll974l: Kaiser! Il984l; 



Bardeen et aflll986l; ICole fc Kaisej 11989 : iBond et al.l 
White! ll996T ) as dueto the fact that the 



1991: IMo 



first galaxies represented rare peaks in the cosmic den- 
sity field. 

As an analogy, imagine searching on Earth for 
mountain peaks above 5000 meters. The 200 such 
peaks are not at all distributed uniformly but instead 
are found in a few distinct clusters on top of large 
mountain ranges. Similarly, in order to find the early 
galaxies, one must first locate a region with a large- 
scale density enhancement, and then galaxies will be 
found there in abundance. The density of stars should 
thus show stron gly biased (i.e., amplifie d) fluctuations 
on large scales (|Barkana fc L ocb 2004). These large- 
scale fluctuations, and their effect on everything from 
feedback to observational predictions, had been pre- 
viously underestimated, in part because the limited 
range of scales available to numerical simulations put 
this topic mostly out of their reach. 

This idea first made a major impact on stud- 
ies o f cosmic reionization. In Bark ana fc Loebl 
(2004) we argued that the typical sizes of H II 
bubbles during reionization should be around 10 
or 20 Mpc (see Figure [2]), while many numerical 
simulations of reionization at the time followed 
a total box belo w this size. Further analyti- 
cal models (jFurlanetto. Zaldarriaga. fc Hernquistl 
2004 aj) and large-scale num e rical simu l ations 
(ICiardi. Ferrara. fc White! 120031: IZahn et all 120071 : 
iMellema et alj 120061 : ISantos et all 12008ft have indeed 
demonstrated the dominance of large bubbles due to 
large groups of strongly-clustered galaxies. This has 
helped motivate the large number of observational 
efforts currently underway, since large-scale fluctua- 
tions are easier to detect (as they do not require high 
angular resolution). 

The same idea of large-scale fluctuations in galaxy 
numbers soon found a n important app li cation in a 
different regime. In Bark ana fc Loebl (|2005bT ) we 
showed that fluctuations in the galaxy number den- 
sity cause fluctuations even in the intensity of long- 




Figure 2: During reionization, the ionized bub- 
bles created by cluste red groups of galaxies 
(Bar kana fc Loebl l2004f) imprinted a signature 
the power spectrum of 21-cm fluctuations 



(JFurlanetto. Zaldarriaga. fc Hernquistl l2004al). 
The i llustration (top panel, from iBarkanal 
( 2006bl )) shows how regions with large-scale 



overdensities form large concentrations of galaxies 
(dots) whose ionizing photons produce large 
ionized bubbles. At the same time, other large 
regions have a low density of galaxies and are 
still mostly neutral. A similar pattern has been 
confirmed in large-scale numerical simulations 
of reionization (e.g., bottom panel, showing a 
two- dimensional slice from a 150 Mpc simulation 
box (|Mellema et al.ll2006f l). 



journals.cambridge.org/pas 



range radiation, leading to a new source of 21-cm 
fluctuations. Specifically, we considered the back- 
ground of Lya radiation. As mentioned above, 
the spin temperature of hydrogen atoms in the 
IGM can be coupled to the gas temperature indi- 
rectly through the Wouthuysen-Field effect, which 
involves the absorptio n of Lya photons. While it 
was previo usly known (Madau, Mciksin, & Roes 1997; 
IChen fc Miralda-Escudei 120041 ) that this Lyq coupling 
likely occurred in the IGM due to Lya photons emitted 
by early stars at z ~ 20 — 30, this radiation background 
had been assumed to be uniform. This intuition was 
based on the fact that each atom sees Lya radiation 
from sources as far away as ~ 300 Mpc. However, 
we showed that relatively large, potentially observable, 
21-cm fluctuations are generated during the era of ini- 
tial Lya coupling, for two reasons: galaxy fluctuations 
are significant even on scales of order 100 Mpc, and 
also a significant fraction of the Lya flux received by 
each atom comes from relatively nearby sources. Since 
a relatively small number of galaxies contributed to 
the flux seen at any given point, we found that Pois- 
son fluctuations could be significant as well, producing 
correlated 21-cm fluctuations. If observed, we argued, 
this new predicted signal would not only constitute the 
first detection of these early galaxies, but the shape 
and amplitude of the resulting 21-cm power spectrum 
would also probe their average properties. 

This has led to a variety of follow-up work, 
including more precise analyses of the a tomic 
cascades of Lyman series photons (|Hiratal 2006; 
IPritchard fc Furlanettd l2006h and a calculation of 
heating due to an in homogcncous X-rays background 
(jPritchard fc Furlanettd 2007) that followed a simi- 
lar approach and found somewhat larger resulting 21- 
cm fluctuations. We have also predicted a significant 
boost in the 21-cm power spec trum from Lya fluc- 
tuations (Naoz & Barkana 200J3) due to the repeated 
scattering of the photons from stars on their way to the 
hydro gen a toms (see also Chuzhov. Alvarez, fc Shapird 
(2006J) and lSemelin. Combes, fc BaeU (|2007l H. 



3 The new game in town: rel- 
ative motion between the 
baryons and dark matter 

As noted above, current observational efforts in 21- 
cm cosmology (and high-redshift astronomy more gen- 
erally) are focused on the reionization era (redshift 
z ~ 10), with earlier times considered much more dif- 
ficult to observe. This has begun to change as a re- 
sult of recent new predictions that suggest that the 
pre-reionization, z ~ 20 era of even earlier galaxies 
is primed for observational exploration. This is made 
possible by a recently noticed effect on early galaxy 
formation that had been previously neglected. 

Up until recently, studies of early structure forma- 
tion were based on initial conditions from linear per- 
turbation theory. However, th ere is an important effect 
missi ng from this treatment (JTseliakhovich fc Hiratal 
2010). At early times, the electrons in the ionized gas 
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Figure 3: The contribution of various scales to 
the mean squared velocity difference between the 
baryons and dark matter (at the same position) at 
recombination. The contribution per log fc of fluc- 
tuations at wavenumbe r k is shown vs. fc. From 
iTseliakhovich fc Hiratal (|2010l) . 



scattered strongly with the then-energetic CMB pho- 
tons, so that the baryons moved together with the pho- 
tons in a strongly-coupled fluid. On the other hand, 
the motion of the dark matter was determined by grav- 
ity, as it did not otherwise interact with the photons. 
Thus, the initial inhomogeneities in the universe led 
to the gas and dark matter having different velocities. 
When the gas recombined at z ~ 1100, it was mov- 
ing relative to the dark matter, with a relative veloc- 
ity that varied spatially. The root-mean-square value 
at recombination was ~ 30 km/s, which is supersonic 
(Mach number 5). 

Figure [3] shows the contribution of fluctuations on 
various scales to the variance of the velocity differ- 
ence. This highlights two important properties of this 
relative motion. First, there is no contribution from 
small scales, so that the relative velocity is uniform in 
patches up to a few Mpc in size; the velocity is gener- 
ated by larger-scale modes, up to ~ 200 Mpc in wave- 
length. The uniformity on small scales is critical as it 
allows a separation of scales between the spatial varia- 
tion of the velocity (on large scales) and galaxy forma- 
tion (on small scales). Each individual high-redshift 
mini-galaxy forms out of a small region (~ 20 kpc for 
a 1O 6 M0 halo) that can be accurately approximated 
as having a uniform, local baryonic wind; the rela- 
tive velocity is thus also referred to as the "streaming 
velocity" . The second important feature of Figure [3] 
is the strong baryon acoustic oscillation (BAO) sig- 
nature. Arising from the acoustic oscillations of the 
photon-baryon fluid before recombination, this strong 
BAO signature is a potential observational fingerprint 
of the effect of this relative motion, as is further de- 
tailed below. 

This relative motion effect is not in itself a surprise, 
but before 2010 it had not been noticed that this effect 
was dropped within the standard approach. The stan- 
dard initial conditions for both analytical calculations 
and numerical simulations have been generated based 
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on linear perturbation theory, in which each k mode 
evolves independently. Indeed, the relative velocity 
is negligible if any single scale is considered. How- 
ever, it is important as an effect of large scales (which 
contribute to the velocity difference) on small scales 
(which dominate early galaxy formation). Specifically, 
the relative motion makes it harder for small-scale 
overdensities in the dark matter to gravitationally ac- 
crete the streaming gas. Now, observing such small- 
scale fluctuations directly would require far higher res- 
olution than is currently possible. Nonetheless, this 
effect is immensely important because of the effect on 
star formation. Since stellar radiation strongly affects 
21-cm emission from the surrounding IGM, 21-cm cos- 
mology offers an indirect probe of the relative velocity 
effect. 

The first effect of the streaming velocity on ha- 
los to be analy zed was the suppression of th e abun- 
dance of halos ( Tscliak hovich fe Hiratall20ld ). Since 
the baryons do not follow the dark matter pertur- 
bations as closely as they would without the veloc- 
ity effect, linear fluctuations are suppressed on small 
scales (where the gravitationally-induced velocities are 
comparable to or smaller than the relative velocity). 
According to the standard theoretical models for un- 
derstandi ng the abundance of halos as a functio n 
of mass (|Press fc Schechterl Il97i iBond et al.lll99ll 1. 
this should result in a reduction of the number den- 
sity of high-redshift halos of mass up to ~ 1O 6 M0 
(|Tseliakhovich fc Hiratall2010h . a mass range that is 
expected to include most of the star-forming halos at 
early times. 

iDalal. Pen fc Seliakl (|2010T ) next pointed out a sec- 
ond effect, namely that separately from the effect on 
the number of halos that form, the relative veloc- 
ity also suppresses the gas content of each halo that 
does form. They also claimed that this second effect 
results in 2 mK, large-scale 21-cm fluctuations dur- 
ing Lya coupling, with a power spectrum showing a 
strong BAO signature due to the streaming velocity 
effect. These conclusions were later shown to be quali- 
tatively on the mar k but invalid quantitatively. In par- 
ticular, wc showed (Tscliakh ovich. Barkana. fc Hirata 
2011; Fialkov et al. 2012) that the gas-content effect is 
a minor one on star-forming halos, and is mainly im- 
portant for the lower-mass gas minihalos that do not 
form stars. Also, a 2 mK signal would in any case be 
too weak to be of major interest to the observers (see 
below). 

Meanwhile, many groups began to run small- 
scale numerical simulations that follow indi- 
vidual collapsin g halos subject to the stream- 



ing velocity (|Ma io. Koop mans. fc C iardi 



Stacv. Bro mm. fc Loebl 



2011|; 



Greif et al 



2011 



2011 



O'Learv fc McQuinnl 120121; iMcQuinn fc O'Learvl 
20121 : iNaoz. Yoshida. fe Gnedinll2012T). In partic ular 
the simulat i ons oflStacv. Bromm. fc Loebl (|201lT ) and 
I Greif et al.l (|201lh indicated the presence of a third 
effect, i.e., that the relative velocity substantially 
increases the minimum halo mass for which stars 
can form from gas that cools via molec ular hydrogen 
cooling. We found (JFialkov et al J 12012 ) a fit to these 
simulation results, and thus developed a general 



model that includes the effect of density as well as all 
three effects of the velocity on star formation. 

Figure [4] illustrates some of the results of the nu- 
merical simulation studies of the effect of the stream- 
ing velocity on galaxy formation. As expected, a larger 
velocity suppresses gas accretion more strongly, in par- 
ticular reducing the amount of dense gas at the cen- 
ters of halos. But beyond just this general trend, the 
relative velocity effect gives rise to very interesting dy- 
namics on small scales. It disrupts gas accretion in an 
asymmetric way, so that filaments of accreting gas are 
disrupted more easily if they are perpendicular to the 
local wind direction. In addition, halos that form in 
regions of relatively high velocity develop supersonic 
wakes as they move through the wind. 

4 Detecting the First Stars at 
Redshift 20 

While numerical simulations are the best, most ac- 
curate method for studying the small-scale effects of 
the relative velocity, they are unable to simultaneously 
cover large volumes. Simulations that successfully re- 
solve the tiny mini-galaxies that dominated star for- 
mation at early times are limited to ~ 1 Mpc volumes, 
and cannot explore the large cosmological scales that 
might be accessible to 21-cm observations (which are 
currently limited to low resolution). 

On the other hand, analytical calculations are 
limited to linear (plus sometimes weakly non-linear) 
scales, and thus cannot directly probe the non-linear 
astrophysics of halo and star formation. Even if the 
results of simulations are incorporated within them, 
analytical approaches assume small fluctuations and 
linear bias (i.e., that the distribution of stars is a pro- 
portionally amplified version of that of the underly- 
ing density), assumptions that break down in the cur- 
rent context, where the stellar density varies by orders 
of magnitude on scales of a few Mpc. Even on 100 
Mpc scales, fluctuations in the gas temperature are as 
large as order unity (see below). There are also addi- 
tional non-linear relationships in this problem such as 
the dependence of 21-cm temperature on gas temper- 
ature. Thus, linear, analytical calculations can only 
yield rough estimates, even for the large-scale fluctua- 
tions. 

Thus, the best way to generate observable pre- 
dictions from this era of early galaxies is with 
a hybrid method that combines linear theory on 
large scales with the results of numerical simu- 
lations on small scales. We recently developed 
such a hybrid method and used it to produce the 
first realistic, three-dimensional images of the ex- 
pected large-scale distribution o f the first st a rs an d 
the resulting 21-cm emission (jVisbal et al.l [2012). 
In our approach we built upon previous hybrid 
meth ods used for hig h-redshift galaxy formation 
(iTseliakhovich fc Hiratal l20ld: IDalal. Pen fc Seliakl 
l201fj : iMesinger. Furlanetto. fc Cenl |2011| 1. We first 
used the known statistical properties of the initial 
density and velocity perturbations to generate a re- 
alistic sample universe on large, linear scales. Then, 
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Figure 4: Effect of relative velocity on individual halos, from numerical simulations (including gravity and 
hydrodynamics). The colors indicate the gas density, which ranges from 10 -26 g/cm (blue) to 10~ 23 g/cm 
(red). Two halos are shown at z = 20, with a total halo mass of 2 x 1O 6 M (top) or 8 x 1O 5 M (bottom). 
Panels show the result for gas initially moving to the right with a relative velocity of (left), 1 (middle), 
or 2 (right) in units of the root- mean-square value of the relative velocit y at z = 20. M. indicates the 
corresponding Mach number at z = 20. From lOTearv fc McQuinnl ( 2012h . 
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we calculated the stellar content of each pixel on our 
grid using the over all model that we had developed 
(JFialkov et al.l 12012'! to describe the streaming veloc- 
ity effect on galaxy formation; this includes various 
analytical models as well as fits to the results of the 
small-scale numerical simulations. 

We assumed standard initial perturbations (e.g., 
from a period of inflation), where the density and ve- 
locity components are Gaussian random fields. Veloci- 
ties are coherent on larger scales than the density, due 
to the extra factor of 1/fc in the velocity from the con- 
tinuity equation that relates the two fields. This is 
clearly apparent in the example shown in Figure \S\ of 
a thin slice of our simulated volume. The density field 
fluctuates on relatively small scales, while the velocity 
field shows a larger-scale cosmic web, with coherent 
structure on scales of order 100 Mpc. This means that 
the largest scales, which are easiest to observe, will be 
dominated by the pattern due to the velocity effect. 

The resulting distribution of stellar density at z — 
20 is shown in Figure [5] Note the large biasing (i.e., 
amplification of fluctuations) of the stars: density fluc- 
tuations ranging up to ±50% yield (without including 
the relative velocities) a field of stellar density that 
varies by over a factor of 20 (both smoothed on the 3 
Mpc scale). The velocity effect produces a more promi- 
nent cosmic web on large scales, marked by large coher- 
ent regions that have a low density of stars, separated 
by ribbons or filaments of high star formation. The ef- 
fect is even more striking at higher redshifts (Figure[7J), 
which substantially alters the feedback environment of 
the first stars. The various types of radiation that pro- 
duce feedback spread out to a considerable distance 
from each source, but this distance is typically not as 
large as the span of the velocity- induced features. This 
means that regions of low velocity (and thus high star 
formation) experience radiative feedback substantially 
earlier than regions of high velocity (low star forma- 
tion). Thus, the substantial effect of the velocities on 
early star formation makes early feedback far more in- 
homogeneous than previously thought. 

Observationally, these degree-scale fluctuations 
will affect various cosmic radiation backgrounds, and 
in particular the history of 21-cm emission and ab- 
sorption, which depends on the timing of the three 
radiative transitions mentioned in Section 11.21 spin- 
temperature coupling due to Lya photons, molecular 
hydrogen dissociation due to LW photons, and heat- 
ing due to X-ray photons. Although the timing is 
still significantly uncertain, the 21-cm coupling due to 
Lya radiation is expected to occur rather early, with 
the X-ray heating fluctuations occurring later and pos- 
sibly overlapping with signi ficant small-halo s uppres- 
sion due to LW radiation (|Visbal et al.l 120121 ). It is 
most promising to focus on the fluctuations due to 
X-ray heating around redshift 20, since these tend to 
produce larger 21-cm fluctuations (and at lower red- 
shifts that feature lower foreground noise) than the 
Lya coupling era [as suggested by previous analyses 
of th ese two periods without the relative velocity ef- 
fect (JBarkana fc Loebl l2005bl ; iPritchard fc Furlanettol 
2007)] . Given the expected timing, we can analyze the 
heating era assuming that Lya coupling has already 




Figure 5: The large-scale density and velocity 
fields in an example of a slice from a simulated 
volum e 384 Mpc on a side (based on IVisbal et al.l 
( 2012J ) , but taken from a different box f rom th e one 
shown in the Figures in IVisbal et al.l ( 20121) . i.e., 
for a different set of random initial conditions). 
The thickness of the slice is 3 Mpc (which is also 
the pixel size of our grid). For the density field 
(top panel), we show the fractional perturbation 
relative to the mean, at z = 20; for the velocity 
field (bottom) , we show the magnitude of the rela- 
tive motion in units of the root-mean-square value 
(the map is independent of redshift in these rela- 
tive units). 
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Figure 6: Effect of relative velocity on the number 
density of stars at redshift 20. For the same slice as 
in Figure [5] we compare the previous expectation 
(top panel) , including the effect of density only, to 
the new prediction (bottom), including the effect 
of the same density field plus that of the relative 
velocity. The colors correspond to the logarithm of 
the gas fraction in units of its cosmic mean value 
in each case. 



Figure 7: Effect of relative velocity on the number 
density of stars at redshift 40. For the same slice as 
in Figure [5j we compare the previous expectation 
(top panel) , including the effect of density only, to 
the new prediction (bottom), including the effect 
of the same density field plus that of the relative 
velocity. The colors correspond to the logarithm of 
the gas fraction in units of its cosmic mean value 
in each case. The color scale spans the same range 
as in Figure [B] for easy comparison. 
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saturated. As for the LW flux, we continue to illus- 
trate here the case of negligible LW feedback (as was 
assumed in Figures [6] and [7J| , but below we bracket 
the effect of the LW flux by also considering the oppo- 
site limiting case where the LW transition has already 
saturated. 

Figure [S] shows the gas temperature distribution at 
z — 20, at the midpoint of the heating transition, i.e., 
when the mean gas temperature is equal to that of the 
CMB. Regions where the gas moved rapidly with re- 
spect to the dark matter (dark red regions, bottom 
panel of Figure [5]) produced fewer stars (dark blue 
regions, bottom panel of Figure [6| and thus a lower 
X-ray intensity, leaving large regions with gas that is 
still colder than the CMB by a factor of several (dark 
blue regions, bottom panel of Figure [8}. The spatial 
reach of X-rays results in a gas temperature distribu- 
tion that is smoother than the the distribution of stars, 
and this brings out the effect of large-scale fluctuations 
and thus highlights the contrast between the effect of 
density and velocity fluctuations. 

During the heating transition, the 21-cm bright- 
ness temperature (Figure [9| mainly measures the gas 
temperature, although it is also proportional to the gas 
density (and to the square root of 1 + z). The form of 
the dependence, Tb ex 1 — TcMB/T gss , makes the 21- 
cm intensity more sensitive to cold gas than to hot gas 
(relative to the CMB temperature). Thus, the large 
voids in star formation produced by a high relative ve- 
locity lead to prominent 21-cm absorption (dark blue 
regions, bottom panel of Figure [9| seen on top of the 
pattern from the effect of density fluctuations. These 
deep 21-cm cold spots are the main observable signa- 
ture of the effect of the relative velocity on the first 
stars. Note that the observed wavelength of this radi- 
ation is redshifted by the expansion of the universe to 
4.4 meters (corresponding to a frequency of 68 MHz) . 

While these Figures illustrate the detailed pattern 
of the effect of relative velocity on the 21-cm inten- 
sity distribution, upcoming experiments are expected 
to yield very noisy maps that likely must be analyzed 
statistically. Figure [10] shows the predicted effect on 
a key statistic, the power spe ctrum of the fluctua- 
tions in 21-cm intensity (from Visb al et al.l IJ2012T O. 
The velocities enhance large-scale fluctuations (blue 
solid curve compared with red dotted), leading to a 
flatter power spectrum with prominent baryon acous- 
tic oscillations (reflecting the BAO signature in Fig- 
ure O. The signal is potentially observable with a red- 
shift 20 version of current instruments (green dashed 
curve). If there is strong LW feedback (solid purple 
curve), then the small galaxies that rely on molecular- 
hydrogen cooling are unable to form; the larger galax- 
ies that dominate in that case are almost unaffected by 
the streaming velocity, so the 21-cm power spectrum 
reverts to the density-dominated shape (compare the 
solid purple and red dotted curves), but it becomes 
even higher since more massive galactic halos are even 
more strongly biased (i.e., clustered). 

Thus, regardless of the timing of the LW feedback, 
the 21-cm power spectrum at the peak of the heating 
transition should feature large observable fluctuations 
on large angular scales. Beyond just detection of the 
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Figure 8: Effect of relative velocity on the gas tem- 
perature Tfc at redshift 20. For the same slice as 
in Figure [5] we compare the previous expectation 
(top panel), including the effect of density only, 
to the new prediction (bottom), including the ef- 
fect of density and relative velocity. The colors 
correspond to the logarithm of the gas (kinetic) 
temperature in units of the CMB temperature at 
z = 20. 
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Figure 9: Effect of relative velocity on the redshift 
20 21-cm brightness temperature TJ, (which mea- 
sures the observed intensity of radio waves emitted 
by intergalactic hydrogen atoms at 21 cm). For the 
same slice as in Figure we compare the previ- 
ous expectation (top panel), including the effect 
of density only to the new prediction (bottom), 
including the effect of density and relative veloc- 
ity. The colors correspond to the 21-cm brightness 
temperature in millikelvin units. 



Figure 10: Signature of the relative velocity in the 
21-cm power spectrum, at the peak of the X-ray 
heating transition at z = 20. We consider the pre- 
diction including the relative velocity effect (blue 
solid curve) or with the effect of densities only (red 
dotted curve) , both for the case of a late LW tran- 
sition for which the LW feedback is still negligible 
at z = 20. These predictions are compared to the 
projected 1-er telescope sensitivity (green dashed 
curve) based on 1000-hour observations with an in- 
strument like the Murchison Wide-field Array or 
the Low Frequency Array but designed t o oper- 
ate at 50-100 MHz (|McQuinn et alj|2006ft . where 



we include an estimat ed degrad ation factor due 
to foreground removal (JLiul 12012)1 : this sensitivity 
is defined as the signal that would yield a mea- 
surement with a signal-to-noise ratio of unity in 
each k bin of size Afe = 0.5fc averaged over an 
8 MHz frequency band. Future experiments like 
the Square Kilometer Array should reach a better 
sensitivity by more t han an order of magnitude 
(|McQuinnet al.ll2006l) . To allow for the possibil- 
ity of feedback, we also show the prediction for the 
opposite limit of maximum feedback, i.e., an early 
LW transition that has already saturated (purple 
solid curve) . In this plot we have fixed the heating 
transition at z — 20 fo r easy c omparison among 
the various cases. From lVisbal et all (|2012|) . 
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signal, only a mild additional accuracy is necessary in 
order to determine whether feedback has suppressed 
star formation in the smallest halos. If it has not, 
than the velocity effect produces strong BAOs on top 
of a flattened power spectrum, in particular raising it 
by a factor of 4 on large scales (k = 0.05 Mpc -1 , wave- 
length 130 Mpc, observed angle 2/3 of a degree) where 
experimental sensitivity is optimal. If this character- 
istic shape is observed it would confirm that million 
mass halos dominated galaxy formation at this early 
epoch. 

In summary, the velocity effect makes it signif- 
icantly easier to detect the 21-cm power spectrum 
during the epoch of the first heating of the IGM. It 
also creates a clear signature by flattening the power 
spectrum and increasing the prominence of the BAOs 
(which are more strongly imprinted in the velocity 
than in the density fluctuations). While Figure [TO] 
considers a single redshift, similar observations over 
the full Az ~ 6 redshift range of significant heating 
fluctuations could actually detect the slow advance of 
the LW feedback process, during which the power spec- 
trum is predicted to continuously change shape, gradu- 
ally steepening as the BAO signature weakens towards 
low redshift. 



5 Open Questions 

As summarized above, the recently-noticed velocity 
streaming effect causes a special form of luminosity 
bias at high redshift. Unlike the standard galaxy bias 
relative to the underlying density field, early galaxies 
were subject to a bias that depends on the velocity 
field. Specifically, the velocity difference between the 
baryons and dark matter affected the smallest star- 
forming halos in which the gas cooled via molecular 
hydrogen cooling. This bias of luminosity relative to 
velocity results in large fluctuations in the abundance 
of the first stars. Since the velocity varies on larger 
scales than the density field, and larger scales are easier 
to observe (since lower resolution suffices to detect the 
corresponding 21-cm fluctuations), the velocity effect 
immensely increases the feasibility for upcoming 21- 
cm observations to detect stars from higher redshifts 
(~ 15 — 25) than previously expected. Moreover, such 
a detection would be both easier to recognize and more 
exciting than previously thought. 

These developments should make the whole pre- 
reionization epoch far more attractive and interesting 
to the whole community of theorists, numerical sim- 
ulators, and observers. In fact, this era is rich with 
little-explored astrophysics. As the first stars formed, 
their radiation (plus emission from stellar remnants) 
produced feedback that radically affected not only the 
IGM but also the character of newly-forming stars. As 
shown above, the velocity effect substantially modifies 
the distribution of the first stellar generations, and in 
fact the star-formation history is changed further since 
radiative feedback (which affects all stars beyond the 
very first ones) itself depends on this distribution. 

In the previous section we have mentioned only 
some first explorations of the radiative feedback dur- 



ing this early era. However, in order to more reliably 
guide observers to the expected 21-cm power spectrum 
from this era, the radiative feedback should be ex- 
plored in much greater detail. For example, the 21- 
cm fluctuations generated during the heating transi- 
tion, studied by IPritchard fe Furlanetto ( 2007) with- 
out the streaming velocity and bv lVisbal et al.l (|2012T ) 
with the velocity included, should be studied together 
with inhomogeneous LW feedback. Also, the 21-cm 
fluctuations during Lya coupling originally studied by 
iBarkana fc Loebl <l2005bh should be studied in the pres- 
ence of the streaming velocity and the LW feedback. 
A full study of these potentially observable signals will 
also require a new generation of numerical simulations 
that will explore the small-scale interaction between 
relative velocity and star formation in the presence of 
various feedbacks such as heating and LW feedback. 

At the same time, significant uncertainty is likely 
to always remain in any prediction that involves com- 
plex astrophysics such as star formation and feedback. 
Thus, only observations can ultimately confirm or ex- 
clude the promising theoretical predictions shown in 
the previous section. We therefore expect increased 
observational efforts focused on this early epoch. Such 
observations would push well past the current frontier 
of cosmic reionization [z ~ 10, t ~ 480 Myr) for galaxy 
searches and 21-cm arrays, and thus would represent 
a new frontier for extragalactic astronomy. 
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